The geodesic rule and the spectrum of the vacuum. 
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Abstract 

We analyze the consequences of a recent argument justifying the vahdity of the "geodesic 
rule" which can be used to determine the density of global topological defects. We derive 
a formula that provides a rough estimate of the number of string-like defects formed in 
a phase transition. We apply this formula to vacua which are spheres. We provide some 
reasons for the deviation of our predictions from the corresponding accepted values. 
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The geodesic rule [1] is a way that allows us to develop predicting schemes for the 
density of global topological defects arising in phase transitions [2], [3]. Two approaches 
exist for justifying the geodesic rule. The first, time-honored and widely accepted, re- 
lies on the minimization of the energy of the system [4], [5]. We recently proposed [6] a 
second approach which relies on the Markovian character of the dynamics and on the 
random orientation of the collisions among the expanding causally connected volumes of 
true vacuum. Such collisions take place inside the meta-stable vacuum background. In 
the present paper we follow the second approach, and propose a formula that provides 
a rough estimate for the number of global topological strings formed in a phase transition. 

Topological strings are of great importance methodologically [2], [3] due to their non- 
perturbative nature. What is even more important, they have severe experimental and 
observational implications [2], [3], [7] in models allowing their formation. Because of their 
great significance, a scheme leading to calculations of their density and distribution such 
as the geodesic rule has to be theoretically justified and its predictions have to be thor- 
oughly tested. Sometime after [1], the importance of the role of the temperature quench 
and the effects of the related critical slowing down of the underlying phase transition were 
pointed out in [8]- [11]. The predictions relying on this modification of the geodesic rule, 
the Kibble-Zurek scenario, are generally considered to be in agreement with experimental 
data [12]-[18] and numerical simulations [19]-[23]. Any discrepancies between such exper- 
imental/numerical data and the theoretical predictions are attributed to the details of the 
dynamics of the particular model under study, rather than being considered a shortcoming 
of the general approach. Despite all these successes, there is still some work that needs 
to be done, especially on the analytical side, in finding formulae expressing the density of 
topological defects. The current study, aims in such a direction. 

Let A4 denote the vacuum of the model under study and (f) & M. be a parametrization 
of it. We arbitrarily pick (f) — on M. as the origin of a normal coordinate system on 
which the calculations rely. As was explained in [6] , the stochasticity of the orientation 
of the causally connected volumes during their collisions and the Markovian character 
of their subsequent coalescence can be described by a probability distribution function 
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P{(j), t) which obeys the Fokker-Planck equation 

^ = DV^PiM (1) 

Here D stands for the diffusion constant of this stochastic process on M. [6]. The 
solution of (1), with the initial condition 

P(0,O) = ^(0) (2) 

is the heat kernel K{0, (j),t) : M x M x R+ ^ R. We recall that K{0, (f), t) has the 
eigenfunction expansion 

oo 

K(O,</.,t) = 5]e-^^V;(O)/,(0) (3) 

where * denotes complex conjugation, Aj, j = 0, 1, . . . are the eigenvalues of on M.. 
Here >R denote the normalized eigenfunctions of on M. with respect 

to the hermitian product 

(/.-, fi) = / /;(<!>) M(l>)d<f> = Sji, J, / = 0, 1, . . . (4) 

M 

induced by the Riemannian metric 5' of on the space of the square integrable 

functions L'^{M). We assume, for simphcity, that the vacuum is a manifold with- 
out boundary. Since is an elliptic, positive semi-definite operator on M., then 

> 0, j = 0, 1, Moreover, since is self-adjoint with respect to the inner product 

(-,■), then fMeR, J =0,1,.... 

We also assume, for concreteness, that the phase transition giving rise to the topolog- 
ical strings is of first order and proceeds by bubble nucleation [6] . Given this assumption, 
consider three such bubbles with corresponding values of the order parameter (f>i = 0, 02, 
phis- When these bubbles collide the geodesic rule states that traverses a piecewise 
smooth loop C on M. Without any loss of generality assume C has base point 
01 = 0. The points on M. at which C is non-smooth, with probabihty almost 1, 
are exactly 0, 02, 03- If C is non-contractible then the collision of the three bubbles 
will give rise to a topological defect. The fact that to determine the formation or not of 
a topological string one must follow a piecewise smooth geodesic loop C, shows that 



2 



the quantity most closely related to such a string formation is not the generic heat kernel 
K{0, (f), t) but its diagonal, namely K{0, 0, t). This can be expressed as 



oo 

K{0,0,t) = Y.e-'^'fj{0)fj{0) (5) 

j=0 



Since all points cf) E M. are equally likely to be the base-point of C, we have to use in 
what follows instead of K{0^0,t) its average over A4. Such an average is called the 
partition function Z(^M,g), a-nd is defined as [24] 

1 °° /■ 

^<--' w - g L miAm (6) 

In (6), Volg{A4.) is the volume of M. with respect to the Riemannian measure induced 
by g. Since the eigenfunctions fj{4>), j — 0,1, . . . are orthonormal with respect to (•, •), 
(6) simplifies to 

^ oo 

An explicit calculation of 2'(_A4,g)(t) amounts to determination of the spectrum Spec(A^) := 
{Xj, j — 0,1, . . .} of the Laplace-Beltrami operator on functions of A4. Conversely, 
one can prove [24] that knowing Z(^j^^g)[t) amounts to determining Spec{Ai). Unfortu- 
nately, it is not known how to explicitly calculate Spec{A4) or equivalently Z(^_\4^g){t) in 
general, except for very few cases of manifolds with very high symmetry as well as some 
of their quotients. It is known, for example [24], that the n-sphere §" endowed with its 
round metric g, has 

Spec(§") = {\k-^k{n + k-l), keN} (8) 
with corresponding multiplicities 

rrik ^ ^j-^^—^ {n + 2k-l), k>l and mo = 1 (9) 
[n — 1)1 kl 

Since [24] 
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VolgiSn = (10) 



we find 



11^ 
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Other examples of spaces that are occasionally found as vacua of classical field theories, 
for which the spectrum is explicitly known [24] are the real and complex projective spaces 
M.P"' and CP" respectively, the tori T", the lens spaces L{p,q), and somewhat more 
rarely, the Heisenberg manifolds. 

Major problems exist in the calculation of the spectrum if {Ai,g) happens to have 
non-positive sectional curvature everywhere. Then, according to the Hadamard-Cartan 
theorem [25], the universal cover of Ai is diffeomorphic to IR^imX from which we cannot 
draw enough useful information that can help determine Spec(A^). In most cases of 
physical interest though, we do not have to worry about manifolds of non-positive sec- 
tional curvature. Indeed, most frequently, the classical/non-thermal vacua A4 happen to 
be compact, semi-simple Lie groups G, or their homogeneous spaces endowed with the 
corresponding Killing-Cartan metrics or their images under Riemannian submersions [25]. 
These spaces have non- negative sectional curvature everywhere [26], so there is no need 
to worry about the implications of the Hadamard-Cartan theorem in future considerations. 

In most other cases we have to settle with considerably less; within the scope of the 
heat-kernel methods employed in this work, the short-time behavior of the trace of the 
heat kernel is probably the most useful approximation. As the word states, such an 
asymptotic expansion is employed when t — > 0+ and states that 

oo 

K{<p,<p,t) ~ (47rt)-^^«.(0,0)f (12) 

i=0 

where Ui{(f),(f)) are polynomials involving the Riemann tensor, its covariant derivatives 
and their contractions at e M.. Substituting (12) into (7), we find the corresponding 
asymptotic expansion of 

oo 

ZiM,g){t) ^a,t^^ (13) 

1=0 

where 

a2i = - — / Ui{(j),(J))^/detg d(j) ieN (14) 
(47r)^~ J 

and 

a2i+i = 0, ieN (15) 
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In particular, 



do = -—a^ Vol,{M) (16) 
(47r)^~ 

«2 = ^ , , \,^M I R(0)Vd^# (17) 
6(47r) 2 7 
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{2||Riem(0)||2 - 2||Ric(<^)f + 5R2((/))} V^d^ ci0 (18) 



where 1 1 • 1 1 denotes the point- wise norm 

||Riem((/.)||2 = Rmp^R"'"', \\RH<I>W = ^p?^"" (19) 

of the Riemann and the Ricci tensors respectively and R denotes the Ricci scalar. The 
coefficient ag is also known [24], but its expression is too long to state and not very 
enlightening. The zeroth order term has a clear geometric interpretation as a multiple of 
Volg{M). The second-order term is proportional to the Euler characteristic x(7Vl) when 
Ai is a Riemann surface. Unfortunately, in all other cases the non-trivial higher order 
coefficients do not have a straightforward geometric or topological interpretation. 
So despite of their usefulness in this asymptotic expansion, they are not of any great 
help in improving our intuition about the behavior of the diagonal of the heat kernel. 
Subsequently their potential physical significance is also relatively limited, so we will not 
expand further upon them. 

During consecutive bubble collisions, piece-wise smooth geodesies of all possible lengths 
on Ai describe the evolution of (p. A more appropriate expression related to the density 
of topological defects is then 



-oo 



2(M,9) = J Z(^M,g){t) dt (20) 


Substituting (7) into (20), we find 



{+00 
f °° 1 





(21) 



In (21), the first term gives an infinite contribution since for a closed manifold like 
the lowest eigenvalue of the Laplace-Beltrami operator is zero Aq = 0. An infinite 
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subtraction of Z{^M,g) is required to make sense out of (21). The second term of (21) can 
be recognized as the value at t = 1 of the zeta function Cv(^) associated with V^, so 
we can write 

The function Cv(^) is meromorphic over C with isolated simple poles at '^""^"■^ G M 
with J = 0, 1, 2, . . . and corresponding residues 

where T{z), 2; G C\{0, — 1, — 2, . . .} denotes the Euler gamma function. If dimM. is 
large enough so that jo — dimM. — 2 e N namely if dimM. > 3 then t — 1 is a 
simple pole of Cv(^) with residue adimM-2- Incase dimM is an odd number, then the 
residue at t — 1 is zero in accordance with the general result stated above. If a vacuum 
for which dimM > 4 happens to also satisfy adimM-2 — —1, and if in addition we 
interpret the improper integral in (22) as the weak limit (Aq = 0) 

+00 +00 

e-^o*(it = lini J e-'^dt (24) 


then the two infinities of (22) cancel exactly each other and the final result for Z{^M.g) is 
finite. In a generic case, however, an infinite subtraction from Z(^j^^g^ is required. For 
the case of spheres S", (11) 



r(^) \^{n + k-2)\ {n + 2k-l) 



■'(S",ff) 



_2 

71+1 



2^-*^ ] g {n-l)\k\ • k{n + k-l) + i [ ^^^^ 




We can check that the integrated partition functions of (25) are, with few exceptions, 
divergent even if we disregard the term corresponding to the zero eigenvalue Aq = 0. To 
proceed further, an infinite subtraction is required. It is not really difficult to understand 
the origin of this infinity. The calculation of Z(^j^^g^ takes into account the contributions 
of geodesic loops of any length. Each geometrically distinct geodesic loop C is counted 
infinite times, since even upon traversing it as many times as we want, it still remains 
a geodesic loop (although non- minimal) . The measure with which each C is weighted 
in Zi^j^^g) is not sufficient to outweigh the contributions of all its multiples, so the total 
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contribution of C to is infinite. The way around such a divergence should now 

be clear: For a general compact A4, its diameter diam{A4) which is defined as 



diam{M) 



sup{d{p,q), p,qeM} 



(26) 



is finite, with d{p, q) indicating the distance of p,q E Ai calculated with respect to g. 
As it was pointed out above, the classical/non-thermal vacua of interest have everywhere 
non-negative sectional curvature. If, moreover, this curvature is everywhere bounded from 
below, away from zero, hj 6 > then according to Myers' theorem [25] diamAA. < 
The curvature condition of Myers' theorem is satisfied by spheres, by projective spaces 
and by the Aloff-Wallach manifolds. Even if JV[ does not fulfill the requirements of 
Myers' theorem though, there is a weaker estimate for dianiAi. If the Ricci, instead 
of the sectional, curvature is uniformly bounded from below, away from zero, by p > 
then Cheng's theorem [25] states that 



Myers' and Cheng's theorems are useful because they provide upper bounds to a global 
property of the vacuum, namely its diameter, which is used in an essential way in the 
rest of the argument. It is also important, from a practical standpoint, that the condi- 
tions of these theorems can be checked through local computations which can be easily 
performed if one knows the metric of the vacuum g. The length of the longest possible 
piecewisc-smooth minimal geodesic loop C that describes the evolution of the order 
parameter /Goldstonc field (p on Ai during the consecutive bubble collisions should 
depend on diamAi. Let such a dependence be denoted by a{diam Ai). It is very 
difficult to find an exact expression for a for a generic At. For the apparently simpler 
case of C being a closed geodesic, it is not even clear that such an upper bound exists. 
However piecewise smooth geodesic loops are less rigid objects than closed geodesies, ex- 
actly because of the almost certain lack of the differentiability of the former at the points 
(pi e Ai, i = 0, 1,2, . . . which are the values that the order parameter attains in the 
colliding bubbles. With such an upper bound on diarn{Ai). only the piecewise smooth 
geodesic loops C with maximum length a^dianiAi) will contribute to the integrated 
partition function Z{^M,g)- In other words, a{diamAi) provides a natural infrared cutoff 
for the allowed lengths of geodesic loops C contributing to Zf^j^^gy On such, "physical 




(27) 
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j Z^M,g){t) dt (28) 



grounds" , (22) reduces to the estimate 

a{diamM) 

Z{M,9) ~ 



Then issues of subtractions of infinity from that arose from the existence of long 

loops on Jv[ no longer persist. Substituting (7) into (28), we get 

In this equation we immediately see that the contribution of the zero cigcnvahic Aq = 
to the sum is a{diamM) ^ which is finite, a fact which justifies some of the comments 
preceding (28). Whether such a scries converges, or is just a formal expression whose 
infinities should be dealt with further, depends on SpecTW. To illustrate this point, 
consider M.i = S^. Recalhng that r(l) = 1 and J2'^=i p ^ T' 

TT 

= Y2 ^^^^ 
By contrast, let A^2 = S^. Then the general term of the series is 

k{k + l) ' 

from which we can show, using the comparison test, that the series diverges. For both of 
the above examples we disregarded the (infinite) contribution of the zero eigenvalue, as 
is usually done in such calculations. 

To determine the density of topological strings, a set of great importance is the set 
of the non-contractible piecewise smooth geodesic loops CM., which is a subset of the 
loop space CM. To get an estimate for the number of string-hke defects, we have to 
calculate the contributions of all such non-contractible piecewise-smooth geodesic loops 
and compare these to Zf^j^^g), namely we want to determine 

dstr = 1^ — / Z,^M,g){t) dt (32) 

Z{M,g) J 

CM 

The expression (32) provides an estimate rather than an exact prediction for the number 
density of topological strings produced in a phase transition. So, we do really expect 
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deviations between the result of (32) and the experimentally measured or numerically 
computed value of the density of topological strings. To highlight this point we notice that 
for A4 = E)^', n > 1 the density is dgtr = as expected due to that 7ri(S") = 0, n > 1. 
For = we have 



This prediction is four times smaller than dgtr — 0.25 predicted for kinks [27], [28] (in 
2 spatial dimensions), and almost thirteen times smaller than dgtr — 0.88 predicted for 
strings in 3 spatial dimensions [19]. Such discrepancies are due to the following: the do- 
main of integration of (32) should be the set of all piecewise smooth loops rather than the 
set of geodesic loops of A^. In reaching (33) we integrated only over the closed geodesic 
loops. We missed therefore a large contribution to dstr coming from the piecewise-smooth 
geodesic loops. In general, the number of piecewise smooth loops over which we have to 
integrate in (32) is far higher than the number of geodesic loops, for any vacuum M.. 
Also, such number of piecewise smooth geodesic loops should increases at a far higher rate 
than the number of geodesic loops, as the length of C increases. Therefore dgtr as was 
calculated in (32) provides a rough estimate about the density of topological strings rather 
than an actual prediction which can be directly compared with experimental data. Such 
numbers will be reasonably close for simple vacua M. with relatively small diameters and 
simple topology. As the metric and topological properties of M. get more complicated 
however, we would expect increasingly larger deviations between the predictions of (32) 
and the experimental data. 

It is probably worth mentioning, at this point, that the dimension of space-time on 
which the model with vacuum Ai is defined is incorporated in the determination of Ai 
itself. To find M. one has to carry out integrations whose value, or even finiteness, de- 
pends crucially on the space-time dimension. To carry out such a calculation at non-zero 
temperature one has to compute the effective potential [29] which takes into account the 
contributions of the thermal fluctuations to the value of the classical potential. The sin- 
gularity structure of the effective potential of a model depends strongly on the space-time 
dimension [29], a dependence which is then inherited to M.. Therefore the dimension of 
the space-time on which the model is defined is used in our considerations only very indi- 
rectly. Certainly none of our arguments rely on it. We naively expect that an increase of 




(33) 



k=l 
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the spatial dimension on which the model is defined will result in large density of defects. 
We can see this, heuristically, by noticing that in higher dimensions there are more bubbles 
available that can collide with any given bubble, so such an increase should be refiected 
by an increase in the density of the formed defects. Whether this argument can be for- 
mally justified or to what extent it is correct could be the subject of another investigation. 
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